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We present a method by which a quantum-mechanical partition function can be approximated 
from below by an effective classical partition function. The associated potential is obtained by a 
simple smearing procedure. For a strongly anharmonic oscillator and a double-well potential, the 
lowest approximation gives a free energy which is accurate to a few percent, even at zero tempera¬ 
ture. 


The path-integral representation of a quantum- 
mechanical partition function 1 
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involves an infinite product of ordinary integrals which, 
after a Fourier decomposition of the periodic paths 
x(t)=x 0 + 2 *=i ( x n e t6>nT +c.c.) with co n = 2rrn/f3 can 
be written as follows: 1 
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The x„ appear implicitly in x (r), the argument of V. If 
we were able to perform the n =£0 integrals, this would 
leave Z as a simple integral 


^ r dx o -fiwu Q ) 

Z ~ J Vw e 


(3) 


In the high-temperature limit we are able to do so. The 
kinetic terms a>\ \x n \ 2 , n=£ 0, develop very sharp minima 
around zero. The x n are always small [of order 
1 /(co n Vj3)] and may be neglected in V(x(r)). The in¬ 
tegrals over dx T f^dx™ become trivial, resulting in the 
well-known classical limit W(x)—> V(x Q ) as T—► oo: 

Z-Z d = J* -^=^exp{ —PV(x 0 )] . 

Because of the similarity between this and (3) we shall 
refer to the function fV(x 0 ) as the effective classical po¬ 
tential 2 and to the integral (3) as the effective classical 
partition function. 


At first sight, introducing JV(x 0 ) may seem of little use 
since the integrals over all x n , n=f= 0, collected in W(x 0 ) 
are, in general, impossible to perform. The purpose of 
this paper, however, is to show that for smooth potentials, 
V(x), there exists a simple way of evaluating these in¬ 
tegrals approximately to a very high accuracy, leading to 
an upper bound for W(x 0 ), to be denoted by W x (x 0 ), and 
to be found according to the following rules. 

(1) Calculate a smeared version of the potential V(x) as 
follows: 
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(2ira 2 ) 1/2 
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(4) 


with an as yet unknown width parameter a 1 . 

(2) Introduce a second parameter ft and form the auxi¬ 
liary potential 


W l (x 0 ,a 2 ,£l)= -^-ln 
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sinh(j8ft/2) 
/9ft/2 


ft 2 2 _ _ 

-—a 2 +V 
2 a 
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(3) Consider a 2 ,ft as functions of x 0 and calculate, at 
each x 0 , the minimum of ^(xo.a^xol.ftlxo)) with 
respect to the parameters a 2 (x 0 ) and ft(x 0 ). The result is 
the desired approximate effective classical potential 


Wi(x 0 )=: min {Jr,(x 0 ,a 2 (x 0 ),ft(x 0 ))) • 

a 2 (x 0 ),ft(jc 0 ) 


( 6 ) 

Explicitly, the minimization with respect to fl gives, at 
each x 0j the following relation between H and a 2 : 
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while the minimization in a 2 determines ft 2 as a function The derivation of these rules starts out with a trial par- 
of the smeared potential: tition function Z\ in which the effect of the potential en- 

^ 2 / x ^ 3 9 2 /ov ergy upon the n=£0 components x n is approximated by a 

n 2 U 0 )=2—7F 2 Uo) = —f k 2 u 0 ) • 

Ofl OXq 


( 8 ) 


Gaussian potential as follows: 


Zi=e F '= f <0x(r)exp 
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where ft 2 (x 0 ) is an arbitrary local curvature of the poten¬ 
tial and L x (x 0 ) is a trial potential depending only on the 
average coordinate x 0 . Both functions Q(x 0 ) and L x (x 0 ) 
will be determined by an extremal principle. 

The ansatz (9) has the virtue that all x n with n^= 0 can 
be integrated out leaving only an integral over x 0 


z t =f 


dx o 


po,{x 0 )/2 -pt,(x 0 ) 
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Moreover, it is straightforward to calculate, within the 
trial partition function Zj, the expectation of the differ¬ 
ence between the true and the trial potential 


F(x(r))- 
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where 


(ID 


<0>i=Z] 1 J &x{t)Osx p 


- fdr 
Jo 


x 2(r) + i^£l [x(T) _ Xo]2 \-p Ll ( Xo) 


2 2 

Indeed, if we use the Fourier representation of the paths and for V(x(r)) a Fourier representation in space, writing 

C *3-v< 
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the expectation of V(x{r)) can be brought to the form 
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in which all x nt n=£ 0, can again be integrated out. The 
result is 

<V(x (t))>! 


considered as a particular case of V(x (r)) and its smeared 
version v a 2 {x ^ x o^ becomes simply yfl 2 U 0 )a 2 (*: 0 ). 

Hence, the expectation (11) reads 


r dxp Pil(xo)/2 -pl^xq) 
1 J ^2^ sinh[0fl(x o )/2] e 


'>(X 0 ) (X 0) > 


_j r dXQ 0ft(x o )/2 —pL^x 0 ) 
1 J sinh[y3fl(x 0 )/2] * 


(15) 


where V 2 , ,(x 0 ) is the smeared potential (4) with 


a 2 (x 0 )=-2 


P n = 1 +^ 2 (*o) 


(16) 


This can be summed up to give the expression (7). 

The Gaussian potential yn 2 (x 0 )[x(r)—x 0 ] 2 can be 


X 


V aHxJ X 0)- 


n 2 (x 0 ) 


a 2 (x 0 )—L l (x 0 ) 


(17) 


The unknown functions ft 2 (x 0 )>£i(*o)> are now ^eter- 
mined by using the extremal principle explained in Ref. 1, 
Sec. 10.3. It is based on the well-known inequality for 
convex functionals which states that the true partition 
function Z is bounded from below by 
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Z>Zj exp 


-(/< 


ClHx 0 ) 


n dr V(x(t))— „ 
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[X(T)-X Q ] 2 

(18) 


Using (10) and (17) and performing the variations in 
DU 0 ) and Li(x 0 ) we see that the best lower bound is ob¬ 
tained when the integrand of (17) vanishes identically in 
x Q . This fixes 

Cl 2 (x 0 ) 2 

£ 1 (*0 ) = V aHx 0 ) (X ^ ~ - 2 - ° *° ’ 


Equation (10) has the form (3) with W( x 0 ) being bounded 
from above locally by W x of Eq. (5). The bound is fur¬ 
ther improved by minimizing W\ with respect to H(x 0 ) 
which yields W x (x 0 ) and thus the result stated in the be¬ 
ginning. 

As an example, let us apply the procedure to the anhar- 
monic oscillator with V(x) = \x 2 + jgx*. The smeared 
potential (4) is 


V a 2=i a2 +Tga*+j(l+3ga 2 )x 2 +jgx* , 


from which we find, via (8), fl 2 =l + 3ga 2 +3gx 2 . We 
now solve Eqs. (7) and (8) at each x 0 by iteration, and cal¬ 
culate W x (xq) from Eqs. (5) and (6). The free energies F { 
resulting from the integral (3) are shown in Fig. 1 and 
compared with the exact F ex , the classical F c \ y and an ear- 



FIG. 1. Different approximate free energies F x (ours); F 0 
(earlier Feynmann’s), F d (classical); compared with F ex (exact), 
for the anharmonic oscillator potential V(x) = \x 2 4* j-gx 4 , as a 
function of /3=l/T, with anharmonic couplings g =0,2,4,40. 
Notice that F d lies far below the exact curve, while F 0 lies far 
above it. The new Fj, on the other hand, fits extremely well, up 
to a few percent, for all g and down to zero temperature. For 
g < 4, the difference between F x and is hardly visible. 


her upper bound F 0 , derived via the method of Ref. 1, 
Sec. 10.3, which corresponds to the special nonminimal 
choice for O, D=0, so Eq. (7) gives a 2 =(3/ 12. 

Even for strong anharmonicity (g =40), the result is 
good down to T =0, up to 3% accuracy. The high quali¬ 
ty of the approximation can be understood by observing 
that, in general, the T —*0 limit of F x is 

Fi — £°=min{FF,U 0 )} 

J T —Q 

= min [K 2 +l/8a 2 j (19) 

a 2 

which is the same as the minimal expectation value of the 

A • A 

Schrodinger operator H= — yd + V(x) in a Gaussian 
wave packet ( 2na 2 )~ 1 /4 exp( 4a 2 x 2 ). For the anharmonic 
oscillator, the minimum is reached at 
a 2= [2( 1 -h3ga 2 ) 1/2 ] _1 . A Gaussian wave packet is 
known to give extremely good ground-state energies for 
many smooth symmetric potentials. Table I compares the 
energies obtained in this way with the known ground-state 
energies of the anharmonic oscillator 3 up to g =4000. 

Another example is the double-well potential 
F(x)= — \x 2 +\ xA +\/4g. When solving Eqs. (7) and 
(8), the quantity fl 2 (x 0 ) may become negative for jc 0 «0. 
This, however, presents no problem since there is always a 
solution n 2 (x 0 ) which remains in the interval 
n 2 G(-47T 2 r 2 ,oo) for which a 2 >0 [see Eq. (16)]. The 
various free energies F x are shown in Fig. 2. They are in 
slightly worse agreement with the exact ones than those in 
Fig. 1, although our approximation is quite reliable up to 
£«5. Notice that at T — 0, the quantum fluctuations 
wipe out the double well for large g >0.36 (see Fig. 3). 
For g <0.36, the two off-centered minima survive. They 
become lower than the central minimum for g <0.325. In 
each case, the absolute minimum gives the position of the 
optimal Gaussian wave packet. Table II compares the en¬ 
ergy E° with the true ground-state energies E° x which 
shows that the worst possible (T =0) error of our approx¬ 
imation is < 16% for g«0.4. 

Our method can easily be extended to the n- 
dimensional Schrodinger problem in which x f =(x)/ is an 
rc-component vector and trial frequency [H 2 (xo)]y in (9) 
an nXn matrix. In the special case that V(x i )= V(x 2 ) is 
rotationally symmetric, we may introduce longitudinal 
and transverse parts of (ijj via 

^• = ^i(Xo)^(tf^Oy/Xo+n^(xg)(5 l7 —*0 i x 0j/*o) 


and W x becomes 


W,= ^~ 


ln 


sinh(y9n L /2) 


P 

+ (n-l)ln 


pn L /2 

sinh(/?H r /2) 


PFLt/2 


— T \^L a L ~*r( n — 1 )n r a r ] -f ^ a 2 a 2(xo) > 


(20) 


where the smeared potential is 
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FIG. 2. Same plot as in Fig. 1, but for the double-well poten- 
tial V(x)=-\x 1 +±gx 4 + 1/4*. 



FIG. 3. Effective classical potential Wi{x) at T= 0 (-) 

and oo (—) [where W\{x)= V{x)] for the double-well potential 
V{x ) = — j-x 2 -f tS x 4 + * at various values of g. 


V °l4 ix J (2 iral) l/1 (2ira]') l/lin ~ x) 

X f d n x'exp{—Y(x'—x) i [alx i Xj/K 2 +aT(8 i j—x i Xj/x 2 )](x'—x)j}V(x' 2 ). 


( 21 ) 


In principle, it is possible to generalize our method and 
treat also a few components x n with N>n >0 more accu¬ 
rately. Using an Ansatz similar to (9), but with trial func¬ 
tions n 2 (n 0 ,... ,x N ), L\(xq, ... ,x N ). The additional 
work, however, proliferates rapidly. The improvement 
proceeds from high to low T and the T—»0 limit is the 
same for any finite N. 

Notice that the Ansatz (9) cannot be improved by al¬ 
lowing the trial function ft 2 (jc 0 ) to be a matrix H 2 n '(x 0 ) in 
the space of Fourier components x n [i.e., by using 
)*»*«' instead of n 2 U 0 )jL \ x n\ 2 ] which 
would also lead to an exactly integrable trial partition 


function. After going through the minimization pro¬ 
cedure, we would fall back upon our solution 
tl 2 m '(x 0 )=8 m 'ClHx 0 ). 

Let us further point out that by inserting an external 
source it is possible to calculate approximate correlation 
functions 


<x(r)x( 0)>! = ^ 


* 0 - 


0n 2 u o ) 


+ G(t) 


>,• 


where G(r) is the periodic temperature Green’s function 
which includes the zero-frequency part 


TABLE I. Comparison of exact ground-state energy (Ref. 3) with the limit E°= lim r _ 0 Fi which is equal to the minimum of 
W\(x Q ) at T —0 and can be obtained by minimizing the energy expectation of a Gaussian wave function centered at x 0 =0. We also 
have listed the energy of the first two excited states E l a and E 2 a . Level splitting to the first excited state is given in column 6. We see 
that it is well approximated by the values of ft(0) at T =0, as it should, due to Eq. (22). 


8 

E° 

E° a 

E l a 

El 

El-El 

H(0) 

fl 2 (0) 

2 

0.7017 

0.696176 

2.32441 

4.327 52 

1.62823 

1.627 

0.2991 

4 

0.8125 

0.803 771 

2.73789 

5.17929 

1.93412 

2.000 

0.2500 

40.0 

1.5313 

1.50497 

5.32161 

10.3471 

3.81694 

4.000 

0.1250 

4000.0 

6.8279 

6.69422 

23.9722 

47.0173 

17.2780 

18.190 

0.0275 
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TABLE II. Same comparison as in Table I, but now for the double-well potential (Ref. 4). At small g, the minimum of the effec¬ 
tive classical potential at T =0 lies at x 0 = ±x min ^O. _____ 


8 

E° 

El 

El 

El 

El-El 

El-El 

X min 

fi(Xfnin ) 

tt {X min ) 

0.1976 

0.650 

0.6192 

0.6730 

1.51 

0.0538 

0.891 

1.943 

1.255 

0.397 

0.4 

0.549 

0.4709 

0.76776 

1.63485 

0.29686 

1.164 

0 

0.486 

1.030 

4.0 

0.598 

0.57728042 

2.08305212 

4.253 57128 

1.505 77 

3.676 

0 

1.634 

0.3059 

40.0 

1.409 

1.377 816 85 

4.995 66652 

9.894742 35 

3.617 85 

8.517 

0 

3.829 

0.1306 


. J_ ^ to.T 1 _1 sinh[(jg—r)H/2] 

o 2 n +(l 2 ~ (3C1 2 sinh(0£l/2) ’ 

whose r=0 value is a 2 + 1/0O 2 . This can be continued 
analytically to the real time retarded Green’s function 

G R (t)=0{t)(x(t)x( 0)>j 

= —/0(f)/ } sin(fl(x 0 )f )\ (22) 

\ il(x 0 ) 1 1 


which in the Coulomb limit m—► 0 reduces to — 1/1 x | 
times the error function erf[ | x | /(2a 2 ) 172 ]. In the latter 
case, the limit (19) becomes 


E°= min 

*2 


= mm 


8a 2 


+ VM) 


1 


8a 2 Vlr (2a 2 ) 172 


showing that, for T =0, the quantity fl(x 0 ) at the poten¬ 
tial minimum x min gives an approximation to the energy 
difference between ground and first excited state. In 
Table I we see that for the anharmonic oscillator this ap¬ 
proximation is quite good. For the double-well potential, 
it is only good at large g. At small g, when the central 
barrier is very high, our approximation lacks the ability of 
describing tunneling. Therefore, tlix^) at T=0 is not 
the small level splitting caused by tunneling but the dis¬ 
tance between the lowest and the second pair of almost de¬ 
generate energy levels. (See Table II.) 

Finally, let us remark that the method is useful also for 
some singular potentials as long as the smearing pro¬ 
cedure makes sense. An example is the Yukawa potential 
V(r)= —e~ mr /r where 




2gmV/2 |,|/ (2a 2)l/2 


Xexp 


t 2 + 



The minimum is reached at a^ in = 9rr/32 and gives 
E°=— 4/377. This is only 15% different from the true 
ground-state energy — \ such that, after subtracting the 
continuum contribution, the effective classical free energy 
F x of a Yukawa potential obtained by our method is, at 
any temperature, more accurate than this. 

We do not want to end without mentioning a number of 
papers 5 which all have added improvements to our first 
attempts 1 at calculating an effective classical potential. It 
appears to us that our approximation W\(xq) is both 
more accurate and easier to handle than any of its prede¬ 
cessors [we were able to use a simple home computer for 
doing the iteration of Eqs. (7) and (8) and the numeric in¬ 
tegration (3)]. It will be interesting to extend the method 
to quantum field theories and to models of statistical 
mechanics. 6 
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